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Statistical entropies of a general relativistic ideal gas obeying Maxwell- 
Boltzmann, Bose-Einstein and Fermi-Dirac statistics are calculated in a general 
axisymmetry space-time of arbitrary dimension. Analytical expressions for the 
thermodynamic potentials are presented, and their behaviors in the high or low 
temperature approximation are discussed. The entropy of a quantum field is 
^ ! shown to be proportional to the volume of optical space or that of the dragged 

O ■ optical space only in the high temperature approximation or in the zero mass 

O [ case. In the case of a black hole, the entropy of a quantum field at the Hartle- 

^ ■ Hawking temperature is proportional to the horizon "area" if and only if the 

qh| horizon is located at the light velocity surface. 
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X ■ 1 Introduction 

H : 

>^ , Recently many efforts have been focused on understanding the statistical origin of black 

hole entropy^ and its interpretation^'^. One of them is the so-called "brick- wall model" 
introduced by t' Hooft^. By using this model, he first showed that the leading entropy of a 
quantum gas of scalar particles propagating outside the event horizon of the Schwarzschild 
black hole is proportional to its horizon area but diverges near the horizon. The divergences 
arise from the infinite one-particle number or state density of levels due to the presence of 
arbitrary high modes near the horizon. To remove this divergence, t' Hooft introduced a 
brick wall cutoff which is related to the horizon only. 
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Many works on the application of the brick-wall model to various kinds of black hole 
have been done for scalar fields^"^^. In four dimensional Schwarzschild black hole, besides 
t' Hooft's pioneer work^, the model has also been used to calculate the entropies for an 
idea gas obeying the usual three kinds of statistics^. They showed that the area law for the 
entropy of a quantum field is due to the quantum statistics. Ghosh and Mitra^ applied it to 
the Schwarzschild dilatonic black hole. In the case of a four dimensional rotating black hole, 
Lee and Kim^'^ demonstrated that the leading term of the entropy of a neutral scalar field is 
proportional to the area of the event horizon and diverges as the system approaches to the 
event horizon. Similar conclusion holds true also for the leading entropy of a complex scalar 
field at the Hartle-Hawking temperature of a Kerr-Newman black hole^. The reason of the 
divergences is attributed to the infinite number of state or the infinite volume of the phase 
space near the horizon. In other words, the origin of the divergence is that the density of 
states diverges at the horizon. In a 2 + 1 dimensional black hole, the entropy of a quantum 
scalar field has also been studied in Ref. 10 and 11. 

It is well known that bosonic fields have a special class of supcrradiant mode solution^^. 
Because its complexity, the contribution to the entropy from supcrradiant modes hadn't 
been considered in Rcf. 7-10. In the case of a Kerr black hole, it is claimed in Ref. 13 that 
the negative contribution to the entropy from supcrradiant modes is divergent in the leading 
order. However, Ho and Kang^^ pointed out that its entropy contribution is positive rather 
than negative and the previous error in Ref. 11 comes from the incorrect quantization of 
the supcrradiant modes. 

It is generally believed now that the entropy of a quantum field in a black hole is pro- 
portional to the area of a black hole horizon but diverges due to the presence of the event 
horizon. However, Alwis and Ohta^^ demonstrated that the free energy and the entropy for 
a scalar field or a Dirac field in the high temperature approximation is proportional to the 
volume of optical space in a static space-time background. Thus the horizon area of a black 
hole must have a certain relation to the volume of its corresponding optical space. 

Then the following questions arise: What relation between the entropy calculated by 
the brick wall model method and that by the optical method? Has any generality among 
the entropies of a quantum field in different dimensional space-times or in different kinds of 
black holes with the same dimensional number? Jing and Yan^^ had studied the entropy 
of a minimally coupled quantum scalar field in the four dimensional general non-extremal 
stationary axisymmetry black hole. However, for a general arbitrary dimensional space-time. 
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there has no similar work presently. To fill up this gap, we investigate quantum statistics of a 
relativistic idea gas obeying three Kinds of the usual statistics: Maxwell- Boltzmann (M-B), 
Bose-Einstein (B-E) or Fermi-Dirac (F-D) in an arbitrary £)-dimensional space-time. 

In this paper, we first calculate the state density of a relativistical idea gas obeying three 
kinds of the usual statistics. We assume that this state density is effective for these three 
kinds of the usual statistics. Then we do thermodynamics in a curved space-time like the 
usual non-relativistical ones^^ in the flat space-time. The paper is organized as follows: In 
Sec. 2, we give a general description of a D-dimensional stationary axisymmetry space-time 
including static space-times and the flat space-time as well as black hole solutions. Then in 
Sec. 3, we derive the constraint on momentum space and evaluate the state density of phase 
space for a given energy. Sec. 4 and 5 is devoted to calculating thermodynamical potentials 
in two cases: the angular velocity of a quantum fleld = and 0, respectively. Ana- 
lytical expressions of the thermodynamical potentials are given, their asymptotical behaviors 
in the high or low temperature are discussed. General discussions about the divergence of 
the entropy of the quantum fields in a black hole background are presented in Sec. 6. Some 
four dimensional space-times arc considered as examples in Sec. 7. Finally, we present our 
conclusions and problems not being considered here. 

2 Description of general space-time 

To begin with, let us firstly give some general description of a D-dimensional stationary 
axisymmetry space-time. As specially important examples, the flat space-time or a static 
space-time as well as black hole solutions are included under our present consideration. 

In general, the line element and electro-magnetic potential of a D-dimensional stationary 
axisymmetry space-time can be expressed in the following form: 



where the metric elements gu, gt<p, gipip and gij are functions of coordinates x — 
{x^,x'^, • • • , x^~^) only. The metric signature is taken as (— , -|-, • • • , -|-) and its submatrix g^j 
is assumed to be diagonal. This space-time (1) has two Killing vectors dt and relating 
to two conserved quantities, energy and azimuthal angular momentum. As examples, it 
contains a large class of black hole solutions and non-black-hole solutions such as the flat 
space-time. Static black holes are included as special case when g^i^ — 0. The event horizon 
f{x) — 0, if it exists, is a null hypersurface determined by 





(2) 
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g''dlf{x) + 2g'^dtd^f{x) + g^%f{x) + g'^d,f{x)d^f{x) = g'^dif{x)d^f{x) = 0. (3) 
Here the contravariant metric elements being 

<7« = g^jV, g'^ = -gtjV, (f^^ = gu/V, g^^ = l/g,^, V = gug^^ - {gt^)\ (4) 

provided the metric determinant is nonsingular, namely, g^ — T> deigij — 'DgD-2 7^ 0. 

The condition that the nontrial null vector N = (iVi, • • • , iVi)_2) exists, where — 
dif{x), is that the sub-determinant (75-2 of the contravariant metric tensor at the horizon 
(r/i) must be equal to zero, namely, (70-2(^/1) — 0- we assert that the metric determinant 
go is non-singular, and so is the contravariant metric determinant gj^^. From the equality 
g]j\2 — T-^dD^j one can know that the location of the horizon is given by the solutions of the 
following combination of equation and inequality: 

9d-2{^) = 0) namely T){x) = 0, and gnix) ^ 0. (5) 

Eq.(2.3) in Ref. 18 is a part of this combination of equation and inequality that deter- 
mines the location of the horizon. For the metric component g^p^p is nonzero at the event 
horizon (r/j), this expression is equivalent to Eq.(4) in Ref. 16, namely, l/gf**(r/i) = 0. 

In general, we can rewrite function T>[x) around the point as 

V{x) = (r - rhYH{x) « (r - r;,)"i/(r/,), (6) 

here r is referred as to a radial coordinate in subspace, H{x) is an analytical function at the 
point Thi a is the order of zeros of function T>{x), or the order of poles of the sub-determinant 
gD-2- In other words, a is the number of duplicate roots of Eq.(5). For the fiat space-time, 
no solution to Eq.(5) exists, so the index a = 0; For a non-extremal black holes, Eq.(5) has 
a single root r/i, so a = 1; In the extremal case, it has a double root r/j, so we have a = 2, 
etc. 

3 Deduction of momentum constraint and state density 

Before calculating the state density of single particle for a given energy from the volume 
of phase space, we firstly derive the constraint on momentum space. We proceed with the 
Lagrange-Hamiltonian method rather than from the semi-classical approximation of Klein- 
Gordon equation. 

The Lagrangian of a relativistic charged particle in the above background space-time (1) 
is given by 
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2>C = gtd'^ + 2gt^i(p + g<p^(p^ + gijx'-x^ + q{Ati + A^ip + Aji'). 



(7) 



Substituting the canonical conjugate momentum given by the following definition: 

Pt = = 9ttt + gt^'^ + qAt = kt + qAt, 

djO. . . . , 
P^' = ^ = Qtft + 9<p<f^ + qA^ ^k^ + qA^, 

Pi^ 9ij^ + = + qAi, 

into the covariant Hamiltonian defined by 7i = ptt + p<^0 + piX^ — C, we get 
2n = - qA,f + 2g'^{p, - qAt){p^ - qA^) + g^^{p^ - qA^f + g'^{p, - qA,){pj - qAj). 

(8) 

Hamilton constraint 2H — —[j^ yields the constraint on momentum space: 

gHI + 2g''^ktk^ + g'^^ + g'^hk^ + /x^ = g'^{p, - qAi){pj - qAj) + /x^ 

+g'\u; + qAtf - 2g'^{u + qAt){m - qA^) + g^'^im - qA^f = 0, (9) 

where we have put pt — —u!,p^ — m, and —kt — cu + qAt. 

Let Pt — dtS,p^ — d^S,Pi — diS{x), where S — —cut + rrup + S{x), we can derive 
Hamilton- Jacobi (H-J) equation which is a rather good semi-classical (W-K-B) approxima- 
tion to Klein-Gordon equation for a complex scalar field ^ = e*"^ in this geometry. 

General speaking, the H-J equation only has well-meaning for a scalar field. However, 
we assume that constraint Eq.(9) works for particles obeying Maxwell- Boltzmann (M-B), 
Bose-Einstein (B-E) or Fermi-Dirac (F-D) statistics, and we will use it to calculate the 
density of single particle in the classical phase space. The reason is that the state density 
evaluated by the classical phase space method is a rather good approximation to degeneracy 
of discrete levels in quantum case, while the latter is, in general, very difficult to be dealt 
with. 

Secondly, we evaluate the density of single particle state in the case that a quantum 
field has a vanishing angular velocity VLq = 0. (For case VLq ^ 0, see below). Momentum 
constraint of Eq.(9) can be recast into form 

^-^ + ^[m-qA^ + ^(a; + qAt)f = ^(c + qAtf - /.^ (10) 
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On the one hand, for a given energy a;, the hypersurface represented by Eq.(lO) is a 
eUipsoid in {D — 1) dimensional momentum space supposed that it satisfies the following 
conditions: 

Qii > 0, — — > 0, > u . 

' --D -gu ^ 

To prevent from appearing an infinite and imaginary state density, these conditions must 
be satisfied. Thus we need to restrict the system in the region that gu < 0. As the metric 
signature is taken as (— , +,•••,+), so we need gij > 0, —gu > 0, then the above conditions 
reduce to 

-,„>0, -V>0, t±^^,^ („) 

-9tt 

The first condition and second one place restriction on the coordinate space, the third 
one gives the lower bound on the range of energy. If these three conditions are satisfied, then 
the volume of phase space is finite. Otherwise, the hyper-surface is noncompact, and the 
state density g{u!) is divergent. The density of states for a given energy is given by taking 
differentiation of phase volume with respect to energy, g{u!) — dr{u!)/{du!), where r(a;) is 
the volume of 2 (D — 1) dimensional phase space for a given energy cu: 



(47r)^r(^)^ ^^-9tt -9tt ^ ^ 

It must be emphasized that a factor 1/2, though not important to the final results, had 
been ignored in many literatures'*'^'^'^^'^^'^^'^"^'^^ in which the authors who used the radial 
wave number to calculate the free energy of a scalar field. The reason is that when one 
takes square roots of the radial wave number, he only takes a positive root and gets rid 
of a negative one. Physically, one discards the negative wave number; Mathematically, he 
gives up another leaf of paraboloid represented by the radial wave number. In the case of 
Minkowski space-time, the phase volume calculated by Eq.(12) is equal to the volume of the 
coordinate space times the volume of the momentum space being divided by a Planckian 
phase factor (2%)^. Thus the total number of single particle state computed from Eq.(12) is 
correct. 

On the other hand, for a given energy u and a fixed azimuthal angular momentum m, 
Eq.(9) represents a compact surface in {D — 2) dimensional momentum space provided it 
satisfies the following angular momentum-energy constraint condition: 
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+ qAt)'' - 2g"^(u; + qAt)(m - qA^) + g^^{m - qA^f + //^ = -g'^hk^ < 0, (13) 

due to Qij > 0, kf > 0. Thus the volume of {2D — 3) dimensional phase space is easily 
computed, 

xl-s''(w + qA,)' + 2j*(w + qA,)(m - qA^) - g'^(m - qA^)^ - (14) 

The state density or the number of modes for a given uj and a fixed m in {2D — 3) 
dimensional space is given by g{uj,m) = dV {uj , m) / {duj) . Using condition (13) and carrying 
out the integration or summation with respect to m, namely, 

r(a;) = — I dmV{uj,m) = V]r(a;,m), 
we can get the same result of the total number of state for a given energy cu 

9{^) = / dmg{u},m) = ^^(u;,m). 
Here and after we take the quantum number m as a continuous variable. 

4 Thermodynamical potential in the case (Qq = 0) 

The volume of phase space in Eq.(12) and its corresponding state density g{u!) are suitable 
to a quantum field that has a vanishing angular velocity but can have a potential $o- It is 
very convenient to use them in the case of static black holes and the flat space-time. We 
assume that a general relativistical bosonic, fermionic idea gas or non-interaction classical 
Boltzmann gas is in thermal equilibrium at temperature in the background space-time 
(1). The free energy for three kinds of the usual statistics is given by 

-EIdu;g{u;,m)e-f^^'^-i^^^), {M - B), 

m 

E / dLug{LU, m) ln(l - e-^('^-'?*o)), {B - E), 

m 

- E / dujg{uj, m) ln(l + e-^('^-«*o)), {F - D). 



(15) 



After carrying the integration by parts on the r.h.s in Eq.(15) and making a substitution 
S — cu + qAt, B — ^0 + Af, the above equation becomes 



-F 



J du}r{uj) < 



g-/3(a;-?$o) 



g/3(w-q*0)+l 



^-P{£-qB)^ {M-B), 

e/3{£-gS)_l ) {B — E), 

ePie-gB)^^, {F - F>). 



(16) 



The substitution £ — —kt — —pt + qAt — uj + qAt corresponds to a gauge transformation 
which doesn't alter the volume of phase space, thus the density of state is invariant under 
such a gauge transformation. So we have r(a;) = r(£^). One can alway choose such a gauge 
$0 that makes B — + At — 0. Substituting the total number of state (namely, phase 
volume) 
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(47r)^r(^) 



~gtt —gtt 



(17) 



into the r.h.s of the second one in Eq.(16), we get the expression for the free energy 

1 ' 



(47r) 



-9d 



dS 



X 



-gtt 



-gtt Jfi\/-gtt 



(B-E) 
{F-D). 



(18) 



If fJ^\/—gtt > qB, then £ > qB, in this case there exists no superradiant mode for a 
quantum bosonic field. However, if ii\/—gtt < qB, then there is an energy interval jJi^f—git ^ 
£ < qB corresponding to the superradiant mode {uj < g$o) as well as another interval 
£ > qB corresponding to the non-superradiant mode {lu > q^o) ■ It is well known that there 
is no superradiant effect for a fermionic field. Because it is somewhat complicated in the 
superradiant case, we here shall not cope with it. Let us suppose £ > iJ>\/—gtt > qB > 0, 
then after some calculation, we get the final results for the free energy in this case. 



V-go 



n=l 



^ e"^^'^(-l)"+^i^D/2(nM/3v^) 
, n=l 



(M-B), 
(B-E), 

(F-D). 



(19) 



Here being the modified Besscl (or MacDonald) function of D/2-th order. 

By using the asymptotic expression of D/2-th order MacDonald function Kd/2{z) at 
small z: 

2Vr(D/2) 



,D/2 







one can obtain the asymptotic behavior of the free energy in the high temperature approx- 
imation (/? — > 0) or in the zero mass case (/x = 0) 
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-F 



E 



n=l 
n=l 



r(D/2)yB- 



D 



ef^qB^ (M - S) 

Cz.(e^^^), (20) 



Here Cois) = E ^ being Riemann zeta function. Vd-i = J d^ '^xd^pJ , ^R^^ is the volume 

n=l " ^ 

of optical space^^ with metric determinant go = QdI {—gtt)^ ■ The metric of optical space is 

df = de + 2^dtd^ + ^ V + ^-^dx'dxK (21) 
9'tt fl'tt 

If we choose such a gauge potential $o that makes B = 0, and introduce a convenient 
statistical factor jV^i = 1, (■£)(!), —(■£)(— 1), for M-B, B-E, F-D statistics, respectively, then 
the free energy will be 

-^-•^O ^z,/2^D ■ (22) 

The free energy in Eq.(22) coincides with that calculated by the heat kernel expansion 
method in Ref. 15. In the case of a Dirac field, they only differs by a factor 2^/^ which is the 
number of a Dirac spinor components in a £)-dimensional space and isn't considered by us 
here. In Ref. 15, the free energy is derived by functional integral method and made a heat 
kernel expansion in the high temperature approximation in the static space-time. However, 
the start-point of ours is a general stationary axisymmetry space-time. 

The entropy is given by 

S - 8^^^ - Un^^^^I^^^ (23) 

The free energy in Eq.(22) and the entropy in Eq.(23) are proportional to the volume of 
optical space and depend on the dimensional number D of the considered space-time. Their 
dependence on the space-time is only related to the temporal component of the metric 
tensor. If the metric tensor gu is nonzero everywhere in the space-time, there exists no 
divergence. However, if g^ vanishes at somewhere, then divergences appear there. Our 
results agree with that in Ref. 15 in the case of a four dimensional space-time. In the 
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case of a black hole, only after introducing a brick wall cutoff and subtracting minor terms, 
can the entropy be proportional to the "area" of the event horizon at the Hartle-Hawking 
temperature 1//?^ = k/{2t^). (This will be illustrated below.) 

5 Thermodynamical potential in the case (Qq 0) 

However, in the case that a quantum field has a nonzero angular velocity VLq in the 
stationary axisymmetry space-time, the matter is slightly different. Zhao and Gui^^ pointed 
out that "physical space" must be dragged by gravitational field with azimuthal angular 
velocity Q.h, and this is also noticed by Lee and Kim in Ref. 7-10 and other author^^. A 
classical relativistic idea gas or a quantum field in thermal equilibrium at temperature 1/(5 
in this background must be dragged too. Therefore, it can be reasonable to assume that the 
quantum field or the classical particle is rotating with an azimuthal angular velocity ^q{x) 
and has a potential $o(^)- For such a modified angular momentum- energy equilibrium 
ensemble^"^ of states of the field, the thermodynamic potential of the system for particles 
with charge q and mass // is given by 



(5W 



-EI du;g{u, m)e-^^"^°-i^°^), {M - B), 

m 

E I dujg{u;,m) \n{l - e-/3('^-'«^o-5*o))^ {B - E), 

m 

- E / dugiu;, m) ln(l + £-^('^-"^^^0-9*0)), (F - D). 



(24) 



Let us define E — qB — u — mJlo — ?^o, B — -\- At -\- O^qA^, then after carrying out 
the integration by parts in the r.h.s of Eq.(24), we obtain 



-W ^ j dET{E) 



^-0(E-qB)^ ^M-B), 

g^(B-9s)_i, (B-E), 



(25) 



el3{E-qB)^l, D). 

The total number of states r{E) is obtained by substituting cu + qAt — E + flQ{m — qA^) 



into Eq.(9) and reducing this equation to Eq.(28) (see below). Now it is expressed as 

1 



r{E) 



(47r)'^ri 



D+l 



~gtt ~9tt 



(26) 



here we have put gu = gtt + '^gt^Qo + g^^o- 

The finiteness of the state density is guaranteed by the following conditions 

E^ 



gu > 0, -V> 0, 



-9tt 



(27) 
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which comes from the compactness of hypersurface determined by 



kikj —gtt 1 



9tip + 9ip<p^o 



9ij -L> 9tt -9tt 



(28) 



To preserve the state density real and finite, we must restrict the system in the region 
that satisfies —gu > because we want —V > as before. This imposes restrictions 



on the angular velocity in the region that Q — y/—V/g^^ < Qq < fl + \/—V/g^^, where 
— —9t(p/9<pif Suppose that E > fJ>y/—gtt > qB > 0, namely, we only consider the case 
that the non-superradiant mode exists for a scalar field, as the calculation is somewhat 
complicated in the superradiant case. Substituting the total number of single particle state 
r{E) into the thermodynamical potential, we have 



W 



1 



(47r)'^r 



-go 



dE 



-Qtt 



E^ 

^-gtt 



] ' 



{B-E) 



(29) 



e/3(is 



The thermodynamical potential for Bose- Einstein statistics coincides with Eq.(19) in 
Ref. 9 for a quantum scalar field in a four dimensional Kerr-Newman black hole geometry. 

Under our assumption that E > ^^/—gtt > qB > 0, we can always select such a gauge 
potential $o that makes B = without violating the above assumption. After carrying out 
the integration with respect to E, we arrive at results 



-W 



I' 



d^-^xdcp 



V-go 



KD/2{fJ'PV-9tt), 



^D/2 



E 

n=l 

^ (-l)"+lKp/2(nM/3v^) 

;;d72 : 

n=l 



(M-B), 

(B-E), 

{F-D). 



(30) 



In the high temperature approximation {(5 — > 0) or in the zero mass case (// = 0), the 
thermodynamical potential has asymptotic behavior: 



W 



T{D/2) 



TT 



D/2 



(/3x/=^)^" 



(31) 



Here Vd-i — I d xdip^-^J^ is the volume of the dragged optical space with determinant 
go — go/i—gtt)^- The metric of the dragged optical space is 
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ds'' = de + 2^^^4^^did(^ + ^ V + ^^dx'dxK (32) 
Qtt 9tt Qtt 

The entropy in the high temperature approximation is given by 

The thermodynamical potential and its corresponding entropy are proportional to the 
volume of the dragged optical space in the high temperature approximation. Apparently, 
they depend on the temperature and the dimensional number of the considered space-time 
as well as the temporal component of the dragged metric. Their divergences count on the 
property of the dragged metric tensor gu- No divergence will appear when gtt is nonzero 
everywhere. When the angular velocity vanishes, they will degenerate to the case con- 
sidered in the last section. 

Using the asymptotic expression of D/2-th order MacDonald function Kd/2{z) at large 
z and taking only its zero order approximation 



oo 



we can obtain the asymptotic behavior of the thermodynamical potential in the low tem- 
perature approximation (/3 — > oo): 



D+l 

n=l n 



(_l)i+ie-"/3M\/-§tt 

2^ 5+1 

_ n=l n 2 

e-f^^^^\ (M - B) 

(B-E) 

-Ci2+i[-e-^'^^^], (F-D) 



0, when f3 ^ oo. (34) 



In the low temperature approximation, the thermodynamical potential and the entropy 

5 — all exponentially tend to become zero suppose that gtt remains finite at every 
point in the space-time. They will be divergent at the point where gu vanishes. This is 
physically reasonable and is consistent with the third law of the usual thermodynamics. 

6 Discussion: horizon and divergence 
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We have stressed that whether the thermodynamical potential W (or the free energy 
F) diverges or not apparently depends upon whether the metric tensor gu (or gu) vanishes 
or not. The divergence appears if and only if gu is equal to zero. Although there exists a 
event horizon on a black hole background, however if gu is nonzero at this horizon, then 
no divergence will appear. Thus a necessary and sufficient condition that the divergence 
appears at the horizon is that the dragged metric tensor gu is equal to zero at this horizon. 
Because gu degenerates to gu when the angular velocity Qq is zero, so we only need to study 
the general case, namely fio 7^ 0. 

Suppose Tc is the p-fold root of equation: gu = 0, then we can recast gu around the point 
Tc into the form 

~9u = (r - r^YGix) ^ (r - r^YGir,), G(r,) = ^-^Jttir,) = ^9^t?\r,), (35) 

pi arP pi 

here G{x) being an analytical function at the point Vc, the {D — 2)-th coordinate x^'"^ = r 
is a "radial" coordinate. In the first one of the above-head equation, the second expression 
is obtained by taking the lowest order approximation. 

Actually the point Vc is located at the light velocity surface^' In the case f^o = 0, 
the surface such that gu = is the infinite red-shift surface. Apparently, the location of the 
light velocity surface depends upon the choice of the angular velocity Qq. Further, let us 
assume^ that the horizon is on the light velocity surface, namely, the location of the horizon 
satisfies equation guij'h) = 0. Near the horizon r/i, the dragged metric tensor gu tends to 
become zero. 

From the asymptotic expression of D/2-th order MacDonald function K£,i2{z) at small 
one can know that the thermodynamical potential near the horizon has the same behavior 
as it does in the high temperature approximation or in the zero mass case: 



- / ^"--^^{^.^o ^ (36) 

Substituting the lowest approximation of the metric tensor gu into the expression of the 
volume of the dragged optical space, the leading behavior of Vd_i near the horizon is given 
by 



In fact, this assumption is reasonable physically in the case of a black hole. 
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[-9ttr 



/ d'^ ^xdip I dr{r — Th) 



^ Dp/2 I -9D 



Here we introduce a small cut-off e and another cut-off L » rh to remove the infra-red 
divergence and the U-V divergence, respectively. The leading behavior of the entropy near 
the horizon is given by 

here ^(r^) = / d^~^xd(fi^ ^Zqd (^ft) is proportional to the "area" of the event horizon. 

The leading entropy of an idea gas obeying three kinds of the usual statistics diverges 
in e^~^f/'^ as the system approaches the horizon of a black hole if and only if the location 
of the horizon is located at the light velocity surface. Under such a circumstance, the 
leading behavior of the entropy at temperature = /t/(27r) is proportional to the horizon 
" area" , but diverges as the brick wall cut-off e goes to zero. The divergence depends on the 
dimensional number D of the space-time as well as the degeneracy p of the horizon. The 
fundamental reason of the divergence is that the volume of the dragged optical space tends 
to become infinite near the horizon which results in that the density of states for a given 
energy E diverges as the system approaches the horizon. 

However, although there doesn't exist a horizon in a non-black-hole space-time, the 
divergence will also appear when the system approaches the light velocity surface. Thus 
our conclusion is that the divergence has no direct relation to the horizon and it is only 
determined by the equation of the light velocity surface. A necessary condition that the 
divergence appears is that the horizon is located at the light velocity surface. In the following 
section, we will use some concrete examples to illustrate this viewpoint. 

7 Examples 

In this section, on the one hand, we will give some concrete examples to discuss the 
divergence problem, on the other hand, we will determine the location of the event horizon 
and its surface gravity in a given black hole geometry. 

Example A: Minkowski space-time {a = 0) 
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In the four dimensional flat geometry, gu — —i,At — 0,A^ — 0, one can select a gauge 
that satisfles B — — 0. If a quantum fleld has a vanishing angular velocity JIq = 0, then 
the total number of states and the free energy are given by 

M^P), (M - B), 

giiiiwffl_ (B-E), (40) 



n=l 



^Urii^, (F-D). 

, n=l 

Here the volume of the optical space V3 = J^^^ ci^a;^/^^ is is equal to the volume of the box 
that confines the idea gas being considered. 

In the low temperature approximation, the free energy and the entropy all tend to become 
zero in e^'^^. In the high temperature approximation, the asymptotic behavior of the entropy 
is given by 

s{n, = Q)^M,^, (41) 

here for convenience, the statistical factor is introduced Af4 = 1, C4(l), — C4(~l)) for M- 
B, B-E, F-D statistics, respectively. The Riemann zeta constants are all known (4(1) ~ 
7r790,C4(-l) = 7/8(4(1) = 77r7720. 

In the spherical coordinates frame, the Minkowski metric is given by ds"^ = —dt^ + dr^ + 
r'^{d9'^ + sin^ Odip^), the metric determinant is = —r^ sin^ 0, the volume of optical space is 

R 

equal to that of the box, namely, V3 — Air J r'^dr — ^R^. 

No horizon exists in the flat space-time for there is no solution satisfying T> — —r"^ sin^ 9 — 
0, g4 — — sin^ ^ 7^ 0, so the indices a — 0, p — 0, the latter due to gu — —1 0. To make 
states density and entropy flnite, one must restrict the size of 3-dimensional sphere, and use 
a box to confine the propagation of a quantum field. The box acts as imposing a boundary 
condition on the quantum field. This will result in quantization of energy and discretization 
of phase space volume as well as state density 

However, when a quantum field has a non-zero angular velocity fio 7^ 0, the dragged 
metric tensor is = — 1 -|- sin^ 9flQ. The thermodynamical potential is given by 
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K,{^if3^,), (M-B), 



oo 



E (42) 

n=l 

g (-l)"+iK2(nM/3A/=g^ ^ (F-D). 
, n=l " 



The entropy at small (5^^J—gtt behaves hke 



4^3 



5(^0 7^0)^^^, (43) 

here the volume of the dragged optical space = 2-n j drdOr^ sin^^(l — sin^ ^fig)"^. 

The entropy diverges as the angular velocity \Qo\ (rsin^)~^. To preserve the volume 
of the dragged optical space and the entropy finite and real, one must restrict the velocity 
|Qo| < (rsin^)~^. Otherwise the entropy will diverge or become imaginary. For a finite fl, 
the hght velocity surface is located at the surface that r sin ^ = ^, then the index p—1. 
Under the condition |Qo'^sin^| < 1, the volume of the dragged optical space is 

rR 



1/3 = 27r V(A; + / r^^+^dr / sin2^+^ OdO 



The condition that the power scries in the above-head expression converges is \^loR\ < 1. 
When the angular velocity Qq = 0, the dragged volume is equal to that of the sphere 

— AttR^/S. When 7^ 0, its zero-th order approximation is also equal to the sphere 
volume. 

Example B: Four dimensional static black hole {a = 1; 2) 

Next let us consider the Reissner-Nordstrom black hole, the metric and electro- magnetical 
potential are 

ds^ = -^df + "^dr^ + r\d9^ + sin^ Odcp^), (45) 
r A 

A = -^dt, A = - 2Mr + (46) 
r 

The metric determinant is 514 = — sin^ 9. Conditions V = —A sin^ 9 = 0, 5^4 = 
r^sin^^ 7^ can be satisfied by the horizon surface equation A = 0. In the non-extremal 
case (M^ 7^ Q^), the metric has two horizons at r± = M ± (M^ — Q^)^^^. Let rh = r±, then 
the horizon function can be rewritten as A = {r — rh)[r — rh + 2{rh — M)] ^ 2{r — rh){rh — M). 
So we have the index a — 1 and function H{rh) — —2{rh — M) sin^ 9. 
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To preserve the state density real and finite, we should restrict the system in the region 
that A > due to —gu > 0. Suppose that the angular velocity is zero, = 0, and choose 
such a potential $o = —-^t — Q/f that makes the chemical potential B = 0, then the light 
velocity surface coincides the the horizon due to gu — — A/r^ = 0. So we have the index 
p = 1 and constant G{rh) = —'2{rh — M)/rl = —2Kh- The leading term of the volume of the 
optical space is given by 

V.3 = 4. / l^^^ ir^—^,J^^ « -J— L_ . _. (47) 

here Ah — 47rr^ is the area of the horizon with the surface gravity being Hh — {vh — M)/r\. 
In terms of the proper distance cut-off from the horizon Vh to + e: 

the volume of the optical space is rewritten as V3 = Ah/ {2n\5'^). The leading behavior of 
the entropy near the horizon is: 

At the Hartlc-Hawking temperature j3 = /t/i/(27r), the entropy of an idea gas is propor- 
tional to the horizon area Ah and diverges in 6~'^: 

5(^lo = 0)«M^. (49) 

In the case of the extremal Reissner- Nordstrom black hole (M^ = Q^), the horizon 
location Vh — Mis the double root of equation A — {r — M)^ — 0, so the indices a — p — 2. 
As the function V ^ -{r - Mfsin^e.gu = -(r - Mf/r'^ k. -(r - MfM''^, so we have 
constant G{rh) = — and function H{rh) — — s\v? 9. The leading term of the volume of 
the optical space near the horizon is given by 

Vs^ Airdr— / dr- — — . (50) 

Near the horizon, the entropy diverges cubically (in e~^): 

5(". = 0)»Ar.3^. (51) 

The entropy of a quantum field in the extremal Reissner-Nordstrom black hole is propor- 
tional to the horizon area Ag — 47rM^ only when the temperature of the system (3 ~ M^/^. 
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The reason why the thermodynamics of the extremal black hole is ill-defined, however, is 
still unclear^^. 

Example C: Four dimensional stationary axisymmetry black hole (a = 1) 

The third example in which we have an interest is the Kcrr-Ncwman black hole. In the 
Boyer-Lindquist coordinates, the metric and the electro-magnetical potential of the Kerr- 
Newman black hole takes the form 

,2 A-a^sin^6l,2 ^r^ + a^-A . o , , 
as — — dt —2 — asm ddtd^p 

+ ^ll±jfl^,^^!^sinHd^^W4 + de% (52) 

Zj za 
A^-^{dt-a sin^ ed(p) (53) 

with the event horizon function A = — 2Mr + + a^, S = + cos^ 9,1) = — A sin^ 9, 
and the metric determinant (74 = — (Ssin^^)^. The location of the horizon = r± — 
M ± (M^ — Q^ — d^Yl"^ satisfies conditions D = 0,5^4 7^ by equation A = 0. For the non- 
extremal case, the index is equal to one (a = 1), and function ii{r\^ — —2{rfi — M) sin^ 9. 

The entropy of a quantum scalar field in the non-extreme Kerr-Newman black hole had 
been considered by many authors^"^'^^ in the case that the scalar field is co-rotating with 
the black hole, namely the angular velocity is a constant Q,h — ^/('^ft + Q^)) ^-^d the potential 

= Qi"h/{i"l + a^), thus the chemical potential tends to become zero near the horizon. 
Other than this choice and the trial choice Qq = 0, however, we choose a local angular 
velocity and a local potential as: 

gt<p _ a{r'^ + - A) 



g^^ (r2 + a2)2 - Aa^sin^^' 



$0 = -{At + nA^) = , ^^^irtx"? -2, - (54) 
^ (r^ + 0^)2 - Aa^ sm^ 9 

The reason why we choose a local velocity and a local potential is that the light velocity 
surface depends upon the choice of the angular velocity. For instance, in the case of — 0, 
the points satisfying — are on the stationary limit surface. The local velocity and 
potential on the horizon are — Qh,^o — ^h, respectively. They tends to become zero at 
infinity. Under such a choice, the chemical potential is always equal to zero {B — 0), and 
the dragged metric tensor becomes: 

gtt = gtt + '2gt^^ + g^^n' = l/y" = -AS/[(r2 + - Aa^ sin^ 9]. 
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The light velocity surface coincides with the horizon under our choice because the equa- 
tion Qu — l/s'** = can be satisfied by A = 0. Thus the index p — 1. To prevent from 
the presence of infinite and imaginary state density, one must restrict the system satisfying 
inequahties: (r^ + a^)^ > Aa^ sin^ ^ > 0. Also, it places a lower bound and an upper bound 



on the angular velocity Hq: f2 — ^J—Vg-^ < < + ^J—Vg-^. If only A > 0, then the 
choice of the velocity f2o = certainly satisfies this restrictions. 

The leading term of the volume of the dragged optical space is given by 

^3 = r dOd^ t drV^.i^r = 271 r de I' drp^^[{r^ + a'f - Aa' sin^ 9]' 

Jo Jrh+e —L> JO Jrh+e LA^ 



Ir^+e —U JO Jrh+e 

sinede [L 7r(r2 + a2)4 Ah 



dr— — ' ~ — ^ X ' ^ ' " arctan 



a? a 



lo rl + a'^ cos^ 9 Jrh+e 2{rh- Myij- - ThY 4efi;^ avh r^' 
The proper distance from to + e is a function of 9: 



rrh+e lY] , f-rh+e 

dr. 



rrh+e / 2j /■'■ 

— \ —dr fti / 

Jrh V A Jn 



rl + a'^cos'^9 2e{rl + cos'^ 9) 



Irn VA Jrh \2{rh-M){r-rh) V - M 

In terms of the proper distance cut-off 5, the dragged optical volume can be rewritten as 

Ah rl + a? cos^ 9 a Ah , . 

~ TTTT? arctan — ^ 77^7^, (a « r/,). 

2k18^ aru rh 2k15^ 

Here Ah — 47r(r^ + a?) is the area of the event horizon, and Kh — {rh — M) / ((r^ + a^) is the 
surface gravity. In the second approximation, we have taken the slow rotating limit. The 
leading behavior of the entropy near the horizon is: 

If we take as the Hartle-Hawking temperature Kh/(2n), the entropy of an idea gas 
near the horizon is proportional to the horizon area Ah and diverges in as S ^ 0: 

This leading behavior of the entropy of quantum fields near the horizon is a general form 
in the black hole background. It is proportional to the horizon area but diverges as the 
system approaches to the horizon. The reason of the divergence is due to the infinite state 
density for a given energy near the horizon. This agrees with the conclusions in Ref. 8 and 
9. 

As a brief summary, we has used the first example to demonstrate that although there 
exists no horizon in the fiat space-time, if a quantum field has a vanishing angular velocity. 
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one must also impose a box on the system to prevent from the presence of infinite state 
density. In the case that a quantum field has a non- vanishing angular velocity, the entropy 
also diverges as the system approaches to the light velocity surface. To guarantee the state 
density finite and real, the angular velocity must be restricted in a certain region. 

The second example has been used to show that the thermodynamics of the extremal 
black hole is very different from that of the non-extremal black hole. In the third one, we have 
chosen a local angular velocity and a local potential other than the popular uniform velocity. 
Both examples have shown that the leading behavior of the entropy of a relativistical idea gas 
near the horizon is proportional to the horizon area and diverges as the system approaches 
the horizon provided that the horizon is on the light velocity surface. 

All examples have illustrated that the entropy of a quantum field is proportional to the 
volume of the optical space or that of the dragged optical space. In the case of a black hole, 
it is proportional to the horizon area only after introducing a brick wall cut-off. In the four 
dimensional black holes, the leading term of the entropy has a common character. 

Other cases can also be considered. In a lower than four dimensional space-time, we have 
assumed that the statistics arc the usual ones. However, this may be problematical. In 2+1 
dimensional planar system, anyons obeys a novel fractional statistics^°, neither the common 
Bose-Einstein statistics nor the well-known Fermi-Dirac statistics. Theoretically, fractional 
quantum Hall effect probably be interpreted by anyonic statistics^^. 

8 Conclusion 

To summarize, a general framework of general relativistical thermodynamics for three 
kinds of the usual statistics has been done in a D-dimcnsional stationary axisymmetry space- 
time. We start from calculating the density of single particle by the classical phase space 
method. The density of single state is invariant under a gauge transformation, however, it is 
suffered by the dragging of the angular velocity. To proceed, we assume that it is effective in 
an arbitrary dimensional space-time for a relativistical idea gas obeying the usual Maxwell- 
Boltzmann, Bose-Einstein or Fermi-Dirac statistics. A particular needed notice is that this 
assumption is probably invalid in a space-time with its dimensional number lower than four. 
In a space-time higher than the usual four dimension, no such problem exists. 

Thermodynamical quantities such as the free energy or the thermodynamical potential 
and the entropy of a quantum field are evaluated. Exact analytical expressions for the free 
energy or the thermodynamical potential are in terms of the modified Bessel functions. In the 
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high temperature approximation, the statistical entropy of a quantum field is proportional to 
the volume of optical space in the case of a vanishing angular velocity of the quantum field. 
This conclusion agrees with the results that the entropy of a quantum field in a static space- 
time background is obtained by heat kernel expansion in the same approximation in Ref. 15. 
In the case that a quantum has a non- vanishing angular velocity, the entropy is proportional 
to the volume of the dragged optical space. In the low temperature approximation, the 
entropy tends to become zero exponentially. 

In general, the entropy of a quantum field in a D-dimensional space-time depends upon 
the temporal component of the metric tensor gu or the dragged metric tensor cju only, as well 
as the dimensional number D. If the dragged metric tensor gu doesn't vanish at every point 
of the space-time being considered, then the entropy is finite provided the system is confined 
by a box. In the case of a black hole, the leading term of the entropy near the horizon is 
proportional to the horizon "area" only when the horizon is located at the light velocity 
surface. The presence of horizon has no direct relation to the divergence of entropy, but it 
introduces an additional brick wall cut-off in place of the restrictions of a box. Although 
there exists a horizon of the black hole, the behavior of the entropy near the horizon will 
also be finite if the light velocity surface doesn't coincides with the horizon. The divergence 
near the horizon has a definite relation to the dimensional number of a space-time. 

As examples, we have discussed the four dimensional entropy of fiat space-time and that 
of black holes. The results agrees with the already-known results in the literatures. Using 
our general formation, one can compare the behaviors among the entropies in different 
dimensional space-times. It might not be a toy of everything, however, we wish it would 
work at least in a higher dimensional space-time. In the case of a bosonic field, we don't 
consider such things as the contribution to the entropy from the superradiant modes and 
a probable existing phenomena of the well-known Bose-Einstein condensation here. In the 
process of our calculation of the free energy or the thermodynamical potential, we have only 
made a small fugacity expansion also, however, we expect to discuss them in other places. 
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